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We consider quantum nonlinear systems with dissipation described within the Caldeira-Leggett 
model, i.e., by a nonlocal action in the path integral for the density matrix. Approximate classical- 
like formulas are derived in order to evaluate thermal averages of observables: the cases of linear 
and nonlinear dissipation are considered, and the framework is extended to the case of many degrees 
of freedom. 



While several applications to condensed matter systems have demonstrated the usefulness of the improved 
effective potential approach, || its suitability to treat open systems was not immediately realized. Previous studies 
were basically aimed to obtain a classical-like expression for the free energy, both in the case of strictly linear Q] , and 
nonlinear dissipation, p 

I \ The effective-potential method, also called pure-quantum self- consistent harmonic approximation (PQSCHA) after 
its generalization to phase-space Hamiltonians , was recently used || to obtain the density matrix of a nonlinear 
system interacting with a dissipation bath through the Caldeira-Leggett (CL) model ||. In the general case the CL 
model starts from the Hamiltonian 
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where (p,q) and (pi,qi) are the momenta and coordinates of the system (here, one degree of freedom) and of an 
| environment (or bath) of harmonic oscillators. The dissipation is said separable if Fi(q) — F(q), while it is said linear 
7—1 . if also F(q) = q. In the path-integral the bath coordinates can be integrated out exactly, leaving the CL euclidean 
7-H ' action 



S[q(u)] = I du r ^q\u)+V{q{u)) - / ^ / du' n(u-u';q(u),q{u')) . (2) 
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For separable and linear dissipation we have 1Z = K(u—u')Q(q(u),q(u')) and 1Z — k(u—u') q(u)q(u'), respectively. 
In the latter case, from Eq. (|l|) one can get the classical equation of motion for the system in the Langevin form, 

Ch . poo 

mq + m / dt' j(t') q(t-t') + V'(q) = , 
O ■ Jo 

. !-h \ and the Laplace transform 'j(z) of the damping function j(t) turns out to be connected with the Matsubara components 
■ of the damping kernel k(u) : 

H ' 

k(u) = m/T 1 J2 n z lUnU k n , K = |i/ n |7(H"n|) > (3) 

where v n = 2nn/ ' (3 are the Matsubara frequencies. We will condider two cases: 
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Ohmic : j(t) — jS(t — 0) , 7(2) = 7 , 

Drude : j(t) = 7 lu d e"" * , 7(2) = 7 w D /(w D + z) , 

where the dissipation strength 7 and the bath bandwidth o>d characterize the environmental coupling; wd — > 00 gives 
the Ohmic (or Markovian) case. 

In the spirit of the PQSCHA, we define a trial action So as the most general quadratic functional with a two-times 
nonlocal term, i.e., 
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So[q(u)] = I du y9 2 ( u ) + w + \m<j? (q(u)-q)' 



— / du k{u-u) (q(u)-q(u')Y 
4 Jo 



(4) 



Here q = (3 1 J q(u) du is the average point of paths and w = w(q), lu 2 — Lo 2 (q), and k(u) — k(u;q), are parameters 
to be optimized by minimizing the right-hand side of the Feynman inequality, F < Fq — T (S — So) s . 

For any observable 0(p,q) the So-average turns out to be expressed || in terms of its Weyl symbol [pi 0(p,q), 
and takes the classical form 



(6) 
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where ((■■■)) is the Gaussian average operating over p and £ with moments 



{(e)) = a(q) = J- ^ 
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0m v 2 + uj 2 {q) + k n (q) fc->o 2mui 

m ^ u 2 (q) + k n (q) 
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the effective potential is defined as V e s(q) = w(q) + a(q) , with 
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The known nondissipative limits are also reported, with / = (3u>/2. 

The variational determination of the parameters gives the following self-consistent equations: 



w(q) = (V(q + 0}-h m ^ 
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du 



^(«;2g + e+(«),e-(«)))}-ifc(«)a-(«)] , (9) 
m W 2 (g) = (V"(q + 0) - 2 jTdu ((^ + %(«;2g + e + («),f_(«)))) , (10) 
mk n (q) = J^du ^(dl + ~dl_)U(u;2q+C+(u),C-(u))^ (1-COS^u) , (11) 
where the last two arguments of 1Z are to be thought as the sum/difference variables q±(u, u') = q(u) ± q(u') and 

(12) 



11 & 1 \\\ - 1 \ 1 ±cosz/„u 

We restrict ourselves from now on to the case of linear dissipation, where k n is given by Eq. (§) and 
«;(?) = (V(q + 0) - \™ 2 a , m W 2 (g) = <V"(g + 0) • 



(13) 



Let us now go over to the application to a many-body system: we consider a y^-chain of particles, described by the 
(undamped) action 



q 1 Pa , fa-g»-i) 2 , n 2 ii 2v 



(14) 
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where a is the chain spacing, fl is the gap of the bare dispersion relation, and g is the quantum coupling. The classical 
continuum model supports kink excitations of characteristic width f2 _1 and static energy ek — 2r2/3g, which is 
naturally used as the energy scale in defining the reduced temperature t = (/3£k) -1 ; an alternative coupling constant 
is defined as the ratio between the quantum harmonic energy f2 and Sk'- Q = 3g/2. We also use the kink length 
in lattice units R = (ail) -1 (R — > oo in the continuum limit). We assume uncorrelated identical CL baths for each 
degree of freedom, and use the low-coupling approximation JtJ in order to deal with the effective potential, by which 
the partition function can be eventually written as 

( 3t 



yAexp j- -^2.[^(qi-qi-i) 2 +v eS { qi ) J, (15) 

«*(«) =ki-u>- ? 2 ) 2 + -d* + § y: ■ (is) 
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The renormalization parameter D = D(t;Q,R;j,u>T>) generalizes Eq. (||) and is the solution of the self-consistent 
equations 

Lul(t) = n 2 [(1 - 3D) + AR 2 sin 2 (fca/2)] , (18) 

where k n = v n j (wc + v n ) and here v n /Q, = 2-nnt/Q. 

We see that while the averages of coordinate-dependent observables tend to the classical behavior (in other words, the 
environment quenches the pure-quantum coordinate fluctuations), those of momentum-dependent ones are enhanced 
due to the momentum exchanges with the environment, resulting in unexpected behaviour for "mixed" quantities like 
the specific heat. 
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FIG. 1. Double-well <p 4 -chain: mean square fluctuation of the site coordinate (qj) vs. t (left), and nonlinear specific heat 
8c(t) — c(t)—c h (t) (right), for Q = 0.2, R = 5. The curves refer to different Drude dissipation strengths 7 (in units of fi), with 
uj-d/Q = 100. Both quantities tend to the classical behaviour for increasing dissipation strength: this is expected for (q~i), while 
it is a non trivial result for 5c(t) = c(t)— Ch(t). 
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FIG. 2. Double- well y> 4 -chain: kinetic (left) and interaction (right) parts of the specific heat, c K (t) and c v (t), for Q — 0.2, 
R = 5. The curves refer to different Drude dissipation strengths 7 (in units of 57), with lj-d/Q. — 100. The kinetic part appears 
to be quenched by dissipation: the kinetic energy (substantially (pi)) becomes indeed larger, but is flatter in t, so its derivative 
c K (i) decreases with 7. 
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